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ABSTRACT 



We present the first spectroscopic study of the globular clusters (GCs) in the 
giant elliptical galaxy (gE) M86 in the Virgo cluster. Using spectra obtained 
in the Multi-Object Spectroscopy (MOS) mode of the Faint Object Camera and 
Spectrograph (FOCAS) on the Subaru telescope, we measure the radial velocities 
for 25 GCs in M86. The mean velocity of the GCs is derived to be = —354^1 
km s _1 , which is different from the velocity of the M86 nucleus (v ga \ = —234 ±41 
km s _1 ). We estimate the velocity dispersion of the GCs, o v = 292I3 2 , km s _1 , 
and find a hint of rotation of the M86 GC system. A comparison of the observed 
velocity dispersion profiles of the GCs and stars with a prediction based on the 
stellar mass profile strongly suggests the existence of an extended dark matter 
halo in M86. We also estimate the metallicities and ages for 16 and 8 GCs, 
>■ I respectively. The metallicities of M86 GCs are in the range —2.0 < [Fe/H] < —0.2 

with a mean value of —1.13 ±0.47. These GCs show a wide age distribution from 
4 to 15 Gyr. 



Subject headings: galaxies: abundances — galaxies: kinematics and dynamics 



OO 

o; 

CN . — galaxies: elliptical and lenticular, cD — galaxies: individual (M86, NGC 4406) 

— galaxies: star clusters 



1. Introduction 

Globular clusters (GCs) are excellent trace rs for studying the formation history of their 



host galaxies ( jLedl2003l ; iBrodie &: Straderl 120061 ). In particular, a giant elliptical galaxy (gE) 
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contains thousands of GCs, located close to the galaxy center to very far away in the outer 
halo. Therefore, GCs can be used as powerful test particles for studying the kinematics and 
chemical evolution of gE halos. 

Up to now, there are seyeral k inemati c studies of the GC systems in nearby gEs: M49 



(|Zepf et al. 



2000 



Cote et all 12003). M60 ([Pierce et al. 



2006 



Bridges et al. 



2006 



Lee et al. 



Cote et al. 



Park et al 



2008 M iHwang et all 120081). M 87 (ICohen & Rvzhov 



200 1 ) . NGC 4636 (ISchuberth et al. 



2010; 



Lee et al. 



Kissler-Patig et al. 



2004a 



1999 



bj; IWoodley et al. 



2010ah. NGC 1399 (Kissler-Patig et al. 



Richtler et al.ll2004j . 120081: 



2006 



1997; IKissler-Patig & Gebhardt 



Chakrabartv fc Ravchaudhur 



20071 ). and NGC 1407 ( Romanowskv et al. 



Schuberth et al 



998 



Minniti et al. 



2010l) . NGC 5128 rtPeng et al. 



1998 
2008 
1998 



20091 ) . The data used for 



these studies were compiled and reanalyzed in iLee et al.l ( )2010al ). who found that the kine- 
matic properties of the GC systems are diverse among the gEs, indicating di verse merging 



and a ccretion histories of gEs (see also the recent study on the M87 GC system (IStrader et al. 
201lh ). 



The re are also several st udies foc using on the spectros copic ages and metallicit ies of the 



gE GCs: 



Cohen et al. 



and 



erai st uaies ioc usmg on tne spectros i 
Jl998l) for M87, beaslev et~al1 J2OO0I) 



(I2010h for NGC 5128. IPierce et all (120061) for M60. 



and 



Kissler-Patig et al.l fll998h and lForbes et al.l J200l|) for NGC 1399. Recently, IPark et al. 



Cohen et al. 



Peng et al. fl2004bh . iBeaslev et al.l d2008h. IWoodley et al.l a nd IWoodley Sz Gomezl 



(bfloj) for M49, 



Cenarro et al.l (120071 ) for NGC 1407, 



( 120121 ) presented a study of spectroscopic ages and metallicities for the GCs in NGC 4636. 
They also compiled the data of all the gE GCs in the literature, and found that the GC 
metallicity distribution in the combined gE sample is bimodal. 

We have been carrying out a project to investigate the spectroscopic properties of GCs 
in nearby galaxies to understand the formation of the GC systems in gala xies. Our s t udy on 
the kinematics of the G C system of the Virgo gE M60 was presented in ILee et al. I J2008bl) 
and Hwang et al. I fhoosi ) and that on the GC system of the spiral galaxy M31 in iKim et al.l 



(£5" and ILee et al.l l2008cl). Recently, we presented measurements of the radial velocities 



fo r the GCs i n NGC 4636 (IPark et al.ll2010l ). and a detailed kinematic analysis of these data 



111 



Lee et al.l (l2010af). We also i nvestigated the chemical properties of NGC 4636 GCs and 



other gE GCs flPark et alJl2012h 



Here, we present a spectroscopic study of the GCs in M86 (NGC 4406), a gE in the 
Virgo cluster. This galaxy is one of the best targets for the spectroscopic study of the GC 
system because it is located close t o the center of th e Virgo cluster and harbors GCs out to 
large radii from the galaxy center ( ILee et al.ll2010bl ). To date, there has been no published 
spectroscopic study of M86 GCs. 
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In contrast to the absence of spectroscopic studies, there have been sev e ral st udies 
of the photometric properties of M86 GCs. For example, iKundu fc Whitmord (j200lf ) and 
Larsen et al.l (j200lf ) found that the color distribution of M86 GCs is bimodal from an analysis 



of HS T/ WFPC2 images. This bimodality was co nfirmed by HST/ACS data (IPeng et al. 



2006h and ground-based wide-field imaging data jRhode k Zepj booi IParkl boij ). The 
radial number density profile of the M86 GC system is approxi mately well fitted by a de 
Vaucouleurs law and power law fehode fc Zeplbooi IParkl l2012h . Basic information about 
M86 is summarized in Table [TJ 

This paper is composed as follows. Section 2 describes the spectroscopic target selec- 
tion, observation, and data reduction. In §3, we identify genuine M86 GCs and list the 
corresponding photometric and spectroscopic data. We present the kinematic properties of 
the M86 GC system in §4, and the metallicities and ages of the M86 GCs in §5. The primary 
results are summarized in the final section. 



Observation and Data Reduction 



2.1. Spectroscopic Target Selection 

We selected the spectroscopic targets from a photometric sample of the GC candidates 
i n M86 iden tified in Washington CT\ images (15.8' x 15.8') taken at the KPNO 4 m telescope 
(|Park!l2012l ). The GCs in M86 appear as point sources in the ground-based KPNO images, 
and we first selected point sources around M86 with colors of 0.9 < (C — Ti) < 2.1 as GC 
candidates. This (C — Ti) color selection criterion is effective in selecting GC candidates 
in early- type galaxi es: the success r ate of t he photometric search for GCs is about 90% 
in th e case of M87 JCote et al.llioOlh . M49 JCote et alJliooi ). and NGC 4636 JPark et al. 
2010l ). We then selected from the bright sources with magnitudes 19 < Ti < 21.5 as the 
spectroscopic targets. A total of 67 targets were chosen, which included the M86 nucleus 
and two known faint galaxies (NGC 4406B, VCC 0833). We also observed two red, bright 
point sources with (C — Ti) ~ 3 to fill the mask gaps. 



2.2. Observation 



We obtained spectra of the 67 targets from observations in the Multi- Object Spec 
troscopy (MOS) mode of Faint Object Camera and Spectrograph (FOCAS; lKashikawa et al. 
20021 ) at the Subaru 8.2 m telescope on 2002 April 21. We observed two circular masks with 
diameters of 6'. Figure [Upresents a grayscale map of the Ti image of M86 taken at the KPNO 
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4 m telescope, which shows the positions of the spectroscopic targets as well as the observed 
masks. We subtracted the M86 stellar light from the original image using IRAF/ELLIPSE0 
task to show the point sources clearly. 

The observational log is given in Table [2j We used a medium-dispersion blue grism 
(300B) with a dispersion of 1.34 A pixel -1 and the order-cut filter L600 covering 3700 — 6000 
A. Seeing during the observation was 0.6". To make the masks, we obtained i?-band pre- 
images with exposure times of 180 s in the FOCAS camera mode under a seeing condition 
of ~ 1.0" on 2002 March 9. Using these ima ges, we constructe d masks with Mask Design 
Pipeline (MDP), a software utility for MOS ( ISaito et al. 1 120031 ). The slit width along the 
dispersion axis was 0.8", and the resulting spectral resolution was R ~ 500. 

We used three 1200 s exposures in mask-C and one 1800 s exposure in mask-1. We 
also obtained the comparison spectra with Th-Ar lamps before and/or after each exposure. 
We used the FOCAS long-slit spectroscopy mode to calibrate the flux, radial velocity, and 
metallicity. We observed a standard star, BD+33d2642, for the flux calibration and five 
Milky Way (MW) GCs (M5, M13, M92, M107, and NGC 6624) for the velocity and metal- 
licity calibrations. We observed the MW GCs with stepping scan mode to sample an area 
larger than that covered by the slit by moving the slit along the dispersion direction. We 
observed these calibration tar gets during the s ame run and also used their spectra in the 
study of the NGC 4636 GCs fjPark et al.l 120101) . Details of the long-slit mode observation 
are given in Section 2.2 and Table 2 of iPark et al.l ( 120101 ) . 



2.3. Data Reduction 

We first applied basic processing techniques (overscan correction, bias subtraction, and 
cosmic ray rejection) to the CCD images using IRAF tasks: the CCD images were ob- 
tained with a pair of 4K x 2K CCDs. We used the FOC 'AS RED /bigimage task in IDL 



( ISaito et al. 1 120031 ) to produce a large single image from a pair of CCD images and cor- 
rect distortions in the optics. We then clipped the two-dimensional spectrum of each target 
out of the single image, and applied a flat-field correction. The spectrum from each two- 
dimensional image was traced, extracted, and sky-subtracted using the IRAF/APALL task. 
We could not extract the spectra of seven faint targets because of a low signal-to-noise (S/N) 
ratio. Wavelength calibration was performed using the Th-Ar lamp spectra with ~40 useful 
emission lines in the 3800 — 6000 A range. The typical rms error for this calibration is ~ 0.8 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy Inc., under contract to the National Science Foundation. 
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A. The flux of the target spectra was calibrated using the flux standard star. 



Sample flux-calibrated spectra of an M86 GC, a foreground star, and the M86 nucleus 
are shown in Figure [2J We classified the target (ID 81) in panel (b) as a K2III star because 
of the broad (~300A) absorption fe ature around the M gb index. This feature is typically 
seen in the K giant star templates (ISantos et al.l 120021 ) but not in GC templates. Several 
absorption lines typically seen in old stellar systems including the G band, H/3, and Mgb are 
clearly visible in the M86 GC spectrum in panel (a). Absorption features in the spectrum 
of the M86 nucleus are much broader than those in the GC and star because of its large 
velocity dispersion. 



3. Membership Determination 
3.1. Velocity Measurement 



We measure d the radial velocities for the targets using the Fourier cross-correlation task, 
IRAF/FXCOR jTonrv fc David Il97flh . A wavelength range of 4200 - 5400 A was used for 
the cross-correlation because of the low S/N at ~ 4000 A and the strong night sky emission 
line [O I] at 5577 A. During the cross-correlation, we fit the continuum of the spectra using a 
spline function with a 2cr clipping for the low level and a 4cr clipping for the high level. Radial 
velocities were measured for each target using the five MW GC templates; velocities for M86 
GCs derived from the M5, M13, M92, and NGC 6624 templates are consistent within la, 
but those from the M107 template differ by 2a. Therefore, an error-weighted average of the 
first four measurements was taken to give the final radial velocity for each target. The error 
in the measured radial velocity is estimated as < t v >= (£e 4 ~ 2 ) -1 / 2 , where e« is the error in 
each measurement. 



The final number of targets with measured radial velocities is 31. The radial velocities 
for 36 objects among the original 67 could not be determined because of the poor quality of 
the spectra. We determined the radial velocity of the M86 nucleus to be v p = —234 ±41 km 



s 1 , w hich is consistent with the previous measurement v p = —244 ±5 km s 1 (jSmith et al. 



20001 ) . We also measured the radial velocities for two faint galaxies: 



777 ± 38 km s" 1 

for VCC 0833 and v p = 949 ± 23 km s~ x for NGC 4406B. These values are also consistent 



with previous measurements (Sloan Digital Sky Survey and IStrauss et al.lll992l ). Errors in 
our velocity measurements range from 20 to 90 km s _1 with a mean error of 49 ± 16 km s _1 . 
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3.2. Membership Determination 

To identify genuine M86 GCs among the 31 objects with measured velocities, we used 
the (C — T\) colors and the radial velocities. The (C — T\) colors are plotted as a function 
of the radial velocity in Figure [3j note that the M86 nucleus and two faint galaxies are 
not included in the plot. The velocity distribution in panel (b) shows that all objects have 
velocities of -900 to +300 km s _1 except for one object with very a large velocity of v p ~ 2400 
km b _1 . We consider the objects with -900 < v p < 300 km s" 1 and 0.9 < (C - T x ) < 2.1 
to be genuine GCs bound to M86 (the velocity of the M86 nucleus is —234 km s _1 ). The 
(C — Tt ) color r ange is equivalent to a metallicity range of —2.24 < [Fe/H] < 0.33 dex 



(ILee et al. Il2008af ). 



Among the 28 objects (excluding the M86 nucleus and two faint galaxies) two do not 
satisfy the GC selection criteria. One target (ID 226) is classified as a star because its color 
is too red, (C — T x ) > 2.1, even though it satisfies the velocity criterion. The other target (ID 
446) satisfies the color criterion but does not satisfy the velocity criterion: v p = 2434 ± 52 
km s _1 . This object seems to be an intracluster GC in the Virgo cluster. 

Because the systemic velocity of M86 is similar to the radial velocities of the MW stars, 
there could be some MW stars that satisfy our GC selection criteria. To remove these stars, 
we performed a careful visual inspection of the target spectra. We found one object (ID 81) 
with stellar spectral features (see Section I2TB1 and Figure 12(b)) and rejected it from the M86 
GC catalog. Thus, 25 genuine M86 GCs out of 28 targets excluding the M86 nucleus and 
two faint galaxies were identified. We note that there may be one or two more stars among 
the 25 GCs considering the selection efficiency for the GC candidates (see Section |2"TT|) . 

The reliability of our M86 GCs was further ch ecked against the ACS Virgo Cluster Sur- 



vey (ACSVCS) source catalog (j Jordan et al.ll2009i ). This catalog provides the classification 



probability for GC candidates in 100 Virgo early-type galaxies. We found that eight GCs 
(ID = 448, 284, 65, 270, 316, 107, 430, and 324) among the 25 M86 GCs are included in 
the ACSVCS catalog, and all of the eight have a GC probability larger than 90%, which 
confirms our classification. We show these eight GCs as open circles in Figure [31 



3.3. A Catalog of the M86 GCs 

Table [3] lists the photometric and spectroscopic data including the metallicities of 16 
GCs using the BH method and the ages and metallicities of eight GCs using the grid method 
(see Section 5). The first column represents identification number. The second and third 
columns give the right ascension and declination (J2000), respectively. The galactocentric 
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radius and position angle are given in columns 4 and 5, re spectively. The magnitude and 
color information in columns 6 and 7 are from iParkl (120121 ) . The eighth column gives the 
radial velocity and its error measured in this study. The ninth and tenth columns give the 
age and metallicity derived from the grid method. The eleventh column gives the metallicity 
derived from the BH method. The final column indicates the corresponding mask from Table 
|2j The M86 GCs are listed first, followed by foreground stars, a probable intracluster GC, 
and the faint galaxies and M86 nucleus. 



4. Kinematics 

4.1. Mean Velocities and Velocity Dispersions of the M86 GCs 

The M86 GCs in our sample are found in a radial range of 42—446" (i.e., 3.4—36.1 kpc, 
se e Figure ED) . The m ean radial velocity for the 25 GCs determined with the biweight location 



of iBeers et al.l (119901 ) is = —354^1 km s _1 , which is smaller than the radial velocity for 
the M86 nucleus (v ga \ = —234 ± 41 km s _1 ). The radial velocity distribution for the M86 
GCs is sh own in the top pane l of Figure HI the distribution appears to be Gaussian. The I 



statistics (iTeague et al. Ill990l ) gave an / value of 1.022, which is smaller than the critical 
value for rejecting the Gaussian hypothesis at the 90% confidence level, 7q.90 = 1.176. This 
suggests that the velocity distribution of the M86 GCs follows a Gaussian distribution. A 
Gaussian fit yields a peak at v p = —343 km s -1 with a width of <j v = 279 km s _1 . 

The bottom panel of Figure H] shows the radial velocities as a function of the projected 
galactocentric radius R. Mean velocities of the GCs in the inner (42" < R < 240") and outer 
(240" < R < 446") regions are similar to that of the GCs in the entire region. The velocity 
difference between the GC samples and the M86 nucleus is more visible in this panel. 



The velocity dispersion of the M86 GCs determined with the biweight scale of lBeers et al. 
(Il990l ) is cr p = 292^32 km s _1 . The velocity dispersion of the GCs in the inner region 



Op = 292^3g km s 1 , is similar to that in the outer region, a p = 3U±l° km s^ 1 . It is noted 
that the velocity dispersion determined in this study could be contaminated by the inclusion 
of possible MW stars because the systemic velocity of M86, v ga x = —234 km s _1 , is in the 
velocity range of the MW stars. 

The plot of the radial velocities as a function of the (C — 7\) colors, T\ magnitudes, and 
position angles (measured from the north to east) in Figure [5] shows that the mean values 
of the radial velocities do not change with either the magnitude or the color. However, 
the mean values of the radial velocities seem to depend on the position angle, showing a 
minimum value at ~ 60° and ~ 300°. This suggests a rotation of the M86 GC system 
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(discussed in detail in the next section). 



4.2. Rotation of the GC system 

In Figure El we show the spatial distribution of the M86 GCs with measured velocities. 
Although the spatial coverage is neither uniform nor large, the spatial segregation of the 
high- and low-velocity GCs relative to the velocity of the M86 nucleus can be seen, which 
indicates a rotation of the GC system (see also Figure [5] (c)). 

The amplitude and axis of rotation for the M86 GC system was measured with the 
following assumptions: (a) the GC system is spherically symmetric, and (b) the rotation 
axis of the GC system lies in the plane of the sky. If the GCs follow any overall rotation, 
the radial velocities will depend sinusoidally on the azimuthal angles. Thus, we can then 
determin e the amplitude and axis of rotation by fitting the radial velocities (v p ) to the 
function (jCote et al.lboOll . l2003i Iflwang et allbood : iLee et al-llioiOah . 



v p (e) 



J sys 



+ (SIR) sin(9 - G ) 



where QR is the rotation amplitude and v sys is the systemic velocity of the GC system. 

Figure [7] plots the radial velocities of the GCs as a function of position angle and an 
overlay of the best fit rotation curve. We use an error-weighted, nonlinear fit of equation ([I]) 
with f sys fixed to the value of the M86 nucleus velocity. This gives a rotation amplitude of 
flR = 228lgQ km s" 1 , which suggests a rotation of the M86 GC system. However, as this 
result is based on a small sample size, it will need to be examined again with a larger number 
of GCs in future studies. 

The orientation of the rotation axis (9 ) is estimated to be 91°^ 2 i> an d this appears to 
be closer to the photometric major axis of M86 (Q p hot — 120°) than to the minor axis. The 
rotation of the M86 GC system around the major axis is consistent with the result based on 
stellar kinematic s, although the spatial coverage for the stellar kinematics was much smaller 



than this study (IKrajnovic et al 



20111 ) 



We summarize the kinematic results of the M86 GC system in Table HJ the range 
of galactocentric radius of the GCs in arcsec, mean value of the radial distance in arcsec, 
number of GCs, mean projected velocity and velocity dispersion about the mean velocity 
(cr p ), position angle of the rotation axis and rotation amplitude estimated using equation 
(DP), velocity dispersion about the best fit rotation curve (o" pr ), and absolute value of the 
ratio of the rotation amplitude to the velocity dispersion about the best fit rotation curve. 
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The uncertainties of the values represent 68% (la) confidence inte rvals that were de termined 



from the numerical bootstrap procedure following the method of ICote et al.l (j200lf ). 



4.3. Existence of Dark Matter Halo 

Here, we investigate the existence of an extended dark matter halo in M86 by a com- 



(Cote et al. 


2001, 


2003; 


Hwane et al. 


2008; 


Lee et al. 


2010a) 



The stellar mass profile is derived from the surface brightness profile of M86, which is then 
used to compute the VDP. 

Assuming that the M86 GC system is sph erically symmetric in the absence of rota- 
tion, we apply the spherical Jeans equation (e.g.. iBinney fc Tremaine 1119871 ) to a dynamical 
analysis of the GC system. The spherical Jeans equation is 



-7- n cX {r)al{r) + ^Mll n d (r)a 2 (r) = -n d (r) 
ar r 



GM 



tot l 



(2) 



where r is the three-dimensional radial distance from the galactic center, n c \(r) is the three- 
dimensional density profile of the GC system, oy(r) is the radial component of velocity 
dispersion, (3 c \(r) = 1 — oJ(r)/c>f (r) is the velocity anisotropy, G is the gravitational constant, 
and M tot (r) is the total gravitating mass contained within a sphere of radius r . The tangential 
component of velocity dispersion, <Je(r), is equal to the azimuthal component of the velocity 
dispersion, cx^(r), in the absence of rotation. The total mass M tot (r) interior out to any 
radius is the sum of the dark matter mass M^ m (r) and stellar mass M s (r). 

To solve equation (J2J), n c \(r) and M s (r) are obtained from the observation. If we fix 
/3 c i(r) as a constant and assume M^ m (r), we can predict the profile of <J r (r). A comparison of 
this with the observed VDP, will provide information about the existence of a dark matter 
halo. Because the observed VDP for the GCs is projected, we also need to compute the 
projected VDP, a p (R): 



N cl (R) 



n cX a 2 r {r) ( 1 



R 



8 R2 

Pel 



r dr 



R 2 



(3) 



where R is the projected galactocentric distance and the surface density profile, N C \(R), is 
the projection of the three-dimensional density profile n c i(r). 
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4.3.1. Density Profile of the M86 GC System 



The three-dimensional density profiles are derived from t he surface number d ensity of 
the M86 GCs u sing the Navar ro- Frenk- White (NFW) profile (INavarro et al.lll997l) a nd th e 



Dehnen profile ( )Dehnenlll993l ). We adopt the surface density profile gi yen in iParkl ( 2012) 



which was derived using data from the HST/ACS Virgo Cluster Survey ( IJordan et al 



20091 ) 



for GCs at R < 2 arcmin and from the KPNO CT\ data for GCs at R > 2 arcmin. The 
surface number density profile N C \(R) is shown in Figure [HJ We fit the surface number density 
profile to the projections of the NFW profile, n c \(r) = no(r/6) _1 (l + r/b)~ 2 and the Dehnen 
profile, n c \(r) = n (r/a)~ r (l +r/a) r ~ 4 . The profile N C \(R) is derived from an integration of 
the three-dimensional density profile n c \{r) as follows: 



N A {R) 



n c \{r) 



r dr 



R 



R 2 



The results of the fit are summarized as follows: 



n. 



n 



NFW, 
cl I 
Dehnen I 



r) 



cl 



r 



0.052 kpc" 3 (r/14.58 kpc) _1 (l + r/14.58 kpc)~ 2 
0.028 kpc" 3 (r/27.94 kpc)- a93 (l + r/27.94kpc)- 3 - 07 . 



(4) 



(5) 



4.3.2. Stellar Mass Profile 



Figure [9] gi yes the rad ial profile of the M86 surface brightness, derived from the KPNO 
i?-band images JParkll2012h and st ellar mass. This surface brightness profile agrees well with 
the results i n iPeletier et al.l (119901 ) that are based on i?-band photometry and with those in 
Caon et al.l (119901 ) that are converted from 5-band photometry. 



A fit of the surface bri ghtness profile of Park! (120121) to the projection of three -dimensional 



luminosity density profile (ICote et al.ll2003l ; iHwang et al.ll2008t iLee et al.ll2010al ) represented 
by 

(3- 7 )(7-2 7 ) L tot n 



4 T\a 

gives 7 = 1.75, L to t = 1-60 x 10 11 Lr,©, and a 
overlaid in Figure [9] (a). 







r 


hoi 







2( 7 -4) 



(6) 



29.94 kpc. The projected best fit curve is 



Figure |H] (b) shows a three-dimensional stellar mass density profile, p s 



derived with an i?-band mass-to-light ratio of T = 6.5 M Q L H 1 Q (determined in the next 
section). From this profile, we obtain an M86 stellar mass profile represented by 



M s (r) 



Attx 2 p s (r)dx = T / Anx 2 j(x)dx 
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T L 



tot 



(r/a) 



1/2 



1 2(3-7) r 



l + (r/a)V2 



(7-2 7 ) + (r/a) 1 / 2 
1 + (r/a)V2 



(7) 



^.5.5. Extended Dark Matter Halo in M86 

If we adopt the stellar mass profile as the total mass profile (i.e., M to t(?") = M s (r)) 
and p s (r) oc j(r) instead of n c \(r), we can derive the radial component of VDPs for the 
stellar system from the Jeans equation. We also assume i?-band mass-to-light ratios (T ) 
and velocity anisotropics of the stellar halo (/3 s (r)). The projected VDPs are then computed 
from the radial component of VDPs through equation fl3]). 

Figure [10] plots the M86 GC VDP measured in this study and the observed M86 stellar 
VDP given in bender et al. ( 1994 ). The stellar velocity dispersion is almost constant around 
220 km s" 1 in the inner region and is smoothly connected to the GC velocity dispersion at 
R « 3 kpc. We also plot the predicted VDPs derived with T = 5.0 M & L^ Q , s (r) = 0.4 

6.5 MqL^q, (3 s (r) = 0.0 (isotropic), which fit the observed stellar 



(radially biased) and T 
kinematic data at R < 2 kpc well 



For comparison, the predicted VDPs for the M86 GCs 



derived with the same stellar mass profile, the GC number density profiles 



n 



NFW/ 
cl \ 



and 



Dchncn , 
cl 1 



r)) determined in the previous section, and several velocity anisotropics (/3 c i(r) 



n 

+0.99 (radially biased), 0.0 (isotropic), and -99 (tangentially biased)) are also shown. The 
VDPs based on the NFW and Dehnen density profiles are not significantly different from each 
other. Note that none of the models agree with the observed VDPs of the GCs at R > 3 kpc, 
which suggests that the mass-to-light ratio is not constant over the galactocentric radius but 
should increase with radius. This clearly suggests the existence of an extended dark matter 
halo in the outer region of M86. 



5. Metallicities and Ages 

5.1. Metallicity and Age Measurement 

We determined the metallicities and ages of the M86 GCs using two methods: (1) 
the BH method, which determines the metallicity through an empiri c al rel ation between 
absorption line index and metallicity developed by iBrodie fc Huchra I ( 119901 ) . and (2) the 
grid method, which derives the metallicity and age from a comparison of the Lick line index 
of the spectrum with that of a single stellar population (SSP) model. We explain each 
method here. 
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5.1.1. BE Method 



Brodie fc Huchra I (119901 ) and iHuchra et al. I (119961 ) presented linear relations between 
absorption line indices obtained from the integrated spectra and mean metallicities to de- 
rive the metal abundances of old stellar systems. Their method was developed to minimize 
systematic effects such as reddening, individual element abundance anomalies, and instru- 
mental effects. They recommend six indices (G band, MgH, Mg2, Fe5270, CNB, and A) as 
primary calibrators among the 12 line indices for the empirical relations. Here, we use only 
four primary line indices to determine the metallicities of the M86 GCs (G band, MgH, Mg2, 
and Fe5270) because of the low S/N of the spectra in the wavelength range for the CNB and 
A indices. 

Each spectrum is first shifted into the rest frame, and we th en measure the absorptio n 
line indices following the prescription of iBrodie &: Huchra I (ll990f ) and lHuchra et al. I ( 119961 ). 
The measured absorption line indices are calibrated to the BH index system with a zero 
point offset, Index(BH) = Index (Subaru) ± of f set, dete r mined from the indices of five MW 
GCs that are common to this study and IHuchra et al. I (119961 ). The offsets we derived are 
0.014 ±0.015 for G band, -0.013 ±0.009 for MgH, -0.021 ±0.018 for Mg2, and 0.008 ±0.009 
for Fe5270. We determine the metallicity from each index as follows: [Fe/H] Gb£md = 11.415 x 
Gband - 2.455, [Fe/H] MgH = 20.578 x MgH - 1.840, [Fe/H] Mg2 = 9.921 x Mg2 - 2.212, and 
[Fe/H] Fe5270 = 20.367 x Fe5270 — 2.086. Finally, we take an error-weighted average of the 
four measurements as the final metallicity value for each GC. Here, the error is the mean of 
the standard deviation. We are able to determine the metallicity for 16 of 25 GCs based on 
this method. The metallicities for the other nine GCs could not be determined because of 
the low S/N of the spectra. 



5.1.2. Lick Index Grid Method 



Lick absorption line indices are useful for determining the metallicity and age of old 
stellar systems from a comparison of the indices derived from the spectra with the line index 



grids predicted from SS P models JTripicco fc Bell1ll995l: llrager et al. II2OO0I: 



2004 



Puzia et al.l 120051 ; iBeasley et al. 



2008 



Trager et al 



200 



Woodlev et al 



this grid method, we adopt the S SP models g i ven b y iThomas et al. 



follow the technique described in iPuzia et al.l (120051 ) and iPark et al. 



romas et al. 
2010h . 



( 200% 12004 



For 



2005|) and 



2012) 



We calibrate our absorption line indices t o the Lick index system as follows. The 
spectra are smoothed with the Lick resolution (jWorthey fc Ottavianil 119971 ) after shifting 
each spectrum into the rest frame. Lick line indices are then derived from the spectra of 
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the M 86 GCs, following the definitions given in IWortheyl ( 119941 ) and IWorthey fc Ottaviani 
( Il997l ). Here, the line index errors are derived from the photon noise in the spectra before 
the flux calibration. The resulting line indices are then calibrated to the Lick system with 
the zero point offset, Jn<iea;(Lick) = Index (Suba ru) + of f set , deter mine d from the spectra 



of five MW GCs that are common to this study, iTrager et al.l ( 119981 ). and Kuntschner et al. 



(120021 ). The offsets we derived are 0.184 ±0.140 for FT3, 165 ±0.187 for M gb, 0.443 ±0.306 
for Fe5270, and -0.103 ±0.162 for Fe5335 (see Table 1 in lPark et all d2012h for the offsets of 
other indices). We determined the Lick line indices for eight GCs in M86, which are listed 
with errors in Table [5j 



The composite index [MgFe]', defined by [MgFe]' = ^/Mgb x (0.72 x Fe5270 + 0.28 x Fe5335), 
is a good metallicity tracer because of its low sensitivity to [a/Fe]. T he index H/3 is an ag e 



indicator and the least sensitive to [a/Fe] among the Balmer lines (IThomas et al.ll2003l ). 
Thus, we determine the metal licity and age of each GC in the H/3 versus [MgFe]' grids pro- 
vided by IThomas et al.l (120031 ) . Figure [TT1 shows the observational indices of the M86 GCs 
in comparison with the SSP model grids for H/3 versus [MgFe]', which indicate SSP models 
with [a/Fe] = 0.2, [Z/H] = -2.25, -1.35, -0.33, 0.0, 0.35, and 0.67 dex, and ages of 0.4, 
0.6, 0.8, 1, 2, 3, 5, 8, 10, and 15 Gyr. All GCs seem to have ages larger than ~5 Gyr 
and metallicities smaller than the solar abundance. For the GCs inside the envelope of the 
model grid, we take the [Z/H] value and age at the nearest model grid interpolated with 
bins of 0.01 dex for [Z/H] and 0.1 Gyr for age. For the GCs outside the envelope, we take 
the val ues of the nearest e nvelope of the model grid in the direction of the error vector as 
done in iPuzia et al.l ( 120051 ). The two outliers with large difference from the model envelope 
along the H/3 axis might be due to a limit of the low resolution integrated spectroscopy or 
of the model grids because these objects are also often shown in studies of MW GCs and 
M31 GCs derived from high S/N spectra (IPuzia et al.l 120021. 120051 ISchiavon et aJboii l as 



well as other gE GCs ( iCenarro et al.l 12007 



Woodley et al.l 120101 ) . To estimate the errors in 



the age and metallicity, we calculate the ages and metallicities of four data points composed 
of H/3 ± error and [MgFe]' ± error in the grid. The difference between the average of these 
four values and the estimate calculated directly from the index is taken as the final error. 

The metallicities of eight GCs derived from the grid method are compared with those 
from the BH method in Figure H2] (a) . The two measurements are broadly consistent within 
the uncertainty. Here, the total metallicity ([Z/H]) derived from the grid method was con - 
vertedinto [Fe/H](grid) using the relation [Fe/H] = [Z/H] -0.94 [a/Fe] (IThomas et al.lliooih . 
For this conversio n, we adopted [a/Fe] = 0.2, which is the mean [a/Fe] for the GCs in gEs 

r 

( 



Park et al.ll2012l ). We also compared the observational {C — 7i) colors of the M86 GCs 



Parkll2012h with th e model colors derived from the [Z/H] and age using the SSP model in 



Marigo et al.l (120081 ). Among the eight GCs, five GCs agree well within their errors as shown 



-14- 



in Figure [T2l (b). 



5.2. Metallicities and Ages of the M86 GCs 



The metallicities and ages of the M86 GCs are listed in Table El and Figure [13] shows 
(a) [Fe/H] versus R, (b) age versus R, (c) [Fe/H] distribution, and (d) age distribution of 
the M86 GCs. It is not easy to derive any systematic trend in the data because of small 
number statistics. However, several features are noted. First, the metallicities of the 16 GCs 
based on the BH method show a wide range of —2.0 < [Fe/H] < —0.2 with a mean value of 
— 1.13 ± 0.47, which is similar to that based on the the grid method (—0.91 ± 0.44). Second, 
metal-rich GCs ([Fe/H] > —0.9) are found only in the inner region (R < 4'). Third, the age 
of the M86 GCs also shows a wide range fro m 4 to 15 Gyr with a mean of 9.7 ± 4.0 Gyr. 
This is similar to the GCs in other gEs (e.g., IWoodley et al.ll2010l ; iPark et al.ll2012l ). 



Figure shows the metallicity as a function of age for the M86 GCs. This figure 
shows a hint for an age-metallicity relation, meaning that the younger G Cs are more metal- 
rich. This relation ship in gEs has also been seen f or the GCs in M60 (jPierce et al.ll2006l ) 
and in NGC 5128 dWoodlev et al.lboiol ). However, IWoodlev etaD (|201oJ) did not strongly 
conclude its existence due to the large biases in their selected GC sample (e.g., extremely 
bright clusters). A more complete analysis with a more comprehensive data set of gE GCs 
is necessary to draw a strong conclusion about the age-metallicity relation. 



6. Summary 

Using the Subaru spectroscopic data of M86 GCs, we studied the kinematic and chemical 
properties of the M86 GC system. Our main results are summarized as follows. 

1. For the first time, we measured the radial velocities of 31 objects in the M86 field: 25 
M86 GCs, two foreground stars, one probable intracluster GC in the Virgo cluster, two 
faint galaxies, and the M86 nucleus. 

2. The mean velocity of the GCs is Up~ = —3541^ km s _1 , which is different from the 
velocity of the M86 nucleus (v g& \ = —234 ± 41 km s^ 1 ). The velocity dispersion of the 
GCs is a p = 292^32 km s -1 . The M86 GC system shows a hint of rotation. 

3. From a comparison of the VDPs predicted from the stellar mass profile with the ob- 
served VDPs of the stars and GCs, we found evidence for the existence of an extended 
dark matter halo in M86. 
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4. We determined the metallicities for 16 GCs using the BH method, and the ages and 
metallicities for 8 GCs using the grid method. The metallicity of the M86 GCs derived 
from the BH method is in the range —2.0 < [Fe/H] < —0.2 with a mean value of 
— 1.13 ± 0.47. The grid method results in similar [Fe/H] values and a mean age of 9.7 
± 4.0 Gyr. 
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Table 1. Basic Properties of M86 



Galaxy 


M v " v aya b 


Reft' c 


e .r./\. mm 


Distance f 


O star 3 


log(Lx) h 


Ngc 1 Sjv J 




(km s -1 ) 


(kpc) 


(deg) 


(Mpc) 


(km s -1 ) 


(erg s" 1 ) 


blue red 


M86 


-22.7 -234 


15.36 


0.33 29 


16.86 


259±28 


42.00±0.002 


1453 968 3.5±0.5 



" V-band absolute total magni tude: iKormendv et al.l |2009). b Systemic velocity: this study . c Eff ective radius in R-band: 
IParkl d2012f) . d Ellipticity at R eff : IParkl j2012h . e Position a ngle of the minor axi s at R eff : IParkl <2012ft . 1 Distance: iMei et al. I 

(2007 ). 9 Mean v elocity stellar dispersion at _R < 45": jB^nder_et_alJ |[l994]) . h Logarithmic value of X- ray luminosity: 
iBeuing et all jl99Sh . i Numbers of blue GCs and red GCs: iRhode fc Zepl J2004h . j Specific frequency of GCs: iRhode fc Zepl 

||2004|1 . 



Table 2. Observing Log for the Subaru FOCAS/MOS Run 



Mask Name 


R.A. 


Decl. 


N(objects) 


T(exp) 


seeing 


Date 




(J2000) 


(J2000) 




(s) 


(") 


(UT) 


Mask-C 


12:26:11.7 


12:56:17 


34 


3x1200 


0.6 


Apr 21, 2002 


Mask-1 


12:25:52.3 


12:59:39 


33 


1x1800 


0.6 


Apr 21, 2002 



Table 3. Radial Velocity, Age, and Metallicity of Globular Clusters in M86 



ID" 


R.A. 


Decl. 


R 


e 


Ti 


(C - TO 


v p 


Age 


[Z/H] 


[Fe/H] B H 


Mask 




(J2000) 


(J2000) 


(arcscc) 


(dog) 


(mag) 


(mag) 


(km s _1 ) 


(Gyr) 


(dex) 


(dex) 





■ i I - - I 1 1 I.I I ' I II I - I -- 



448 


12 


26 


07 


53 


12 


58 


37 


6 


127.1 


331.0 


21 


31 


± 


0.02 


1 


41 


± 


03 


2 ± 51 
















c 


284 


12 


26 


07 


63 


12 


58 


10 


4 


103.3 


324.4 


20 


88 


± 


0.02 


1 


28 


± 


02 


-441 ± 59 






-1 


20 


± 





29 


c 


65 


12 


26 


13 


87 


12 


57 


57 


7 


77.7 


23.5 


19 


72 


± 


0.01 


1 


49 


± 


02 


-551 ± 24 


7.3 ± 2.4 


-0.89 ± 0.18 


-0 


98 


± 





25 


c 
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26 


05 


54 


12 


57 


36 


8 
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299.0 


20 


84 


± 


0.01 


1 


41 


± 


02 


-46 ± 47 


4.5 ± 3.6 


-0.50 ± 0.34 


-1 


07 


± 





2.3 


c 
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26 


03 


70 
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57 


18 


1 
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20 


98 


± 


0.01 


1 


32 


± 


02 


— 763 ± 64 
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26 


15 


94 


12 


57 


08 





65.0 


70.6 


20 


09 


± 


0.01 


1 


04 


± 


02 


-517 ± 26 


6.6 ± 3.8 


-0.35 ± 0.21 


-0 


44 


± 





26 
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12 


26 


14 


59 


12 


56 


55 


3 
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21 
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± 


0.02 


1 
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± 


05 
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Table 4. Kinematics of the M86 Globular Cluster System 



R (R) 
(arcsec) (arcsec) 


N W p 

(km s _1 ) 


a p ©o 
(km s -1 ) (deg) 


ttR 

(km s _1 ) 


(km s" 1 ) 


42 - 446 203 


25 -3541^ 


9Q9+ 32 91 + 19 


228tgo 


oco+36 n oi+0.32 
Z ° Z -33 U - 8i -0.30 



Table 5. Lick line indices and errors 



ID 


H/3 


Mg2 


Mgb 


Fe5270 


Fe5335 




(A) 


(mag) 


(A) 


(A) 


(A) 


65 


2.423 ± 0.340 


0.120 ± 0.006 


2.765 ± 0.340 


0.996 ± 0.359 


1.074 ± 0.459 


270 


2.527 ± 0.641 


0.075 ± 0.012 


3.316 ± 0.623 


2.496 ± 0.659 


-2.172 ± 0.909 


107 


2.142 ± 0.458 


0.179 ± 0.008 


3.635 ± 0.444 


2.037 ± 0.457 


-0.147 ± 0.602 


143 


0.427 ± 0.474 


0.102 ± 0.008 


3.255 ± 0.448 


2.762 ± 0.478 


2.387 ± 0.574 


58 


2.173 ± 0.301 


0.109 ± 0.005 


0.814 ± 0.307 


2.114 ± 0.312 


2.031 ± 0.399 


79 


2.164 ± 0.368 


0.072 ± 0.007 


0.916 ± 0.363 


2.585 ± 0.380 


5.069 ± 0.475 


96 


2.487 ± 0.312 


-0.005 ± 0.006 


1.128 ± 0.314 


3.116 ± 0.319 


1.920 ± 0.420 


149 


1.705 ± 0.432 


0.043 ± 0.007 


1.164 ± 0.404 


1.955 ± 0.430 


0.134 ± 0.556 
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Fig. 1. — Grayscale map of the T\ image of M86 taken with the KPNO 4 m telescope. The 
M86 stellar light is subtracted from the original image to highlight the point sources. The 
large circles and plus sign indicate the observed masks and the center of M86, respectively. 
The small circles, diamond, triangles, and squares represent the GCs in M86, a probable 
intracluster GC in the Virgo cluster, background galaxies, and foreground stars, respectively. 
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Fig. 2.— Sample spectra : (a) a GC in M86 (ID=284) with T\ = 19.72 mag, (C-7\) = 1.49, 
and [Fe/H] = -0.99 dex, (b) a star (ID=81) with T\ = 19.89 mag classified as a K2III red 
giant star, and (c) the M86 nucleus. All spectra are plotted in the rest frame, smoothed 
using a boxcar filter with 6.7 A, and normalized at 5870 A. 
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Fig. 3. — (C — Ti) color as a function of the radial velocity of the M86 GCs : (a) radial 
velocity distribution, (b) (C — Ti) versus radial velocity and (c) (C — Ti) distribution. 
The circles, crosses, and square represent the M86 GCs, foreground stars, and a probable 
intracluster GC, respectively. The open circles indicate the GCs included in the ACSVCS 
catalog. The box indicated by the dashed line represents the boundary for selecting the M86 
GCs. The vertical dotted line indicates the radial velocity of the M86 nucleus. 
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Fig. 4. — (a) Radial velocity distribution and (b) projected galactocentric radii versus radial 
velocities for the M86 GCs. The large open squares represent the mean radial velocities of 
the GCs in each radial bin represented by the long vertical error bars. The horizontal error 
bars denote the velocity dispersions in each radial bin. The vertical dot-dashed and dotted 
lines indicate the mean velocity and velocity of the M86 nucleus, respectively. The solid 
curved line in (a) is a Gaussian fit of the data. 
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Fig. 5.— Radial velocities as a function of (a) (C — Ti) color, (b) T% magnitude, and 
(c) position angle, 0, for the M86 GCs. The large open squares indicate the mean radial 
velocities of the GCs in each bin, represented by a long horizontal error bar. The vertical 
error bars denote the velocity dispersions of the GCs in the radial bins. The horizontal 
dotted line indicates the radial velocity of the M86 nucleus. 
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Fig. 6. — Spatial distribution of the 25 identified M86 GCs. GCs with velocities larger and 
smaller than the velocity of the M86 nucleus are plotted by open circles and filled boxes, 
respecti vely. The solid ellipse repre sents the boundary for the standard major diameter D 2 § 
of M86 ( jde Vaucouleurs et al.lll99ll ). The dashed lines represent the photometric major and 



minor axes. 
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Fig. 7. — Radial velocity versus position angle of the M86 GCs. The solid curve is a best 
fit rotation curve from Table 4. The horizontal dot-dashed line represents the velocity of the 
M86 nucleus. The vertical arrows show the photometric minor axis of M86. 
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Fig. 8. — Projected surface number density profile for the M86 GCs. The filled and open cir- 
cles represent the G C candidate s from the HST/ACS images and those from t he KPNO CT\ 
i mage s, respectively (jParkll2012l ). The small squares represent the results from lRhode fc Zepi 
( 120041 ) . The solid and dashed lines indicate the projected best fits using the NFW and Dehnen 
density profiles, respectively. 
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Fig. 9. — (a) Surface bri ghtness prof iles of M86. The open and fi lled circles are .R-band sur- 
face photometry given in iParkl (120121 ) and in lPeletier et al.l (119901 ) . respectively. The crosses 
repres ent the surface photometry converted from I?-band photometry given in ICaon et al. 
(11990l ). The dotted line represents a projected best fit of the three-dimensional luminosity 
density profile, (b) Three-dimensional stellar mass density profile derived with a constant 
-R-band mass-to-light ratio of T = 6.5 M Q L^ Q . 
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Fig. 10. — VDPs fo r the stars and GCsin M86. The open squares represent the stellar veloc- 
ity dispersions from bender et all ( 1994 ). The filled circles show the velocity dispersions of 
the GCs. The solid line represents the stellar VDP expected from the stellar mass model with 



a constant stellar mass-to-light ratio of T = 6.5 M L^ Q and a stellar velocity anisotropy 
of (3 S = 0.0. The triple-dot-dashed curve shows the stellar VDP with T = 5.0 M Q L~^~ & and 
P s = 0.4. Also shown are VDPs expected from the same stellar mass model as above but 
with NFW GC density profiles and velocity anisotropics of (3 C \ = 0.0 (short dashed line) and 
+0.99, —99 (long dashed lines). The dotted lines represent the VDPs with a Dehnen density 
profile at /3 d = +0.99 and -99. 
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Fig. 11. — H/3 versus [MgFe]' for the M86 GCs. The grids represent the SSP models with 
[a/Fe] = 0.2 for va rious values of [ Z/Hl (-2.25, -1.35, -0.33, 0, 0.35, and 0.67) and ages (0.4 
-15 Gyr) given by iThomas et al.l (120031 ). 
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Fig. 12. — (a) Comparison of the M86 GC metallicities measured using the BH and grid 
methods, (b) Comparison of the photometric colors ob tained from the observation and a 
model. (C — Ti)o (obs) are the observational colors from iParld (12012 ) and (C — T\)q (model) 
ar e the model c o lors d erived from the [Z/H] and age of the M86 GCs using the SSP model 
of iMarigo et all ( 120081 ). The solid lines represent one-to-one relations. 
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Fig. 13. — Metallicities and ages of the M86 GCs as a function of galactocentric radii, and 
their histograms. The filled circles and hatched histograms indicate the values from the grid 
method, while the open squares and solid histogram are those from the BH method. The 
horizontal dot-dashed and dotted lines indicate the mean values of each parameter measured 
from the grid and BH methods, respectively. 
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Fig. 14. — Age versus metallicity of the GCs in M86. The circles and squares represent the 
GCs of M86 and NGC 5128, respectively. 



